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The existence of a nonperturbative intrinsic heavy 
quark component in the nucleon is a rigorous predic¬ 
tion of Quantum Chromodynamics (QCD). An unam¬ 
biguous experimental confirmation is still missing and 
would represent a major discovery. The goal of this ar¬ 
ticle is to summarize our current understanding of this 
subject with a particular focus on the potential of a high 
energy and high luminosity fixed-target experiment us¬ 
ing the LHC beams (AFTER@LHC) |IH1] to search for 
intrinsic charm. 

Production processes sensitive to the intrinsic heavy 
quark distributions of protons and nuclei are among the 
most interesting hadronic physics topics that can be in¬ 
vestigated with AFTER@LHC. In contrast to the famil¬ 
iar extrinsic contributions which arise from gluon split¬ 
ting in perturbative QCD, the intrinsic heavy quarks 
have multiple connections to the valence quarks of the 
proton and thus are sensitive to its nonperturbative 
structure. For example, if the gluon-gluon scattering box 
diagram, gg —>■ QQ —>■ gg (the analog of QED light-by¬ 
light scattering), is inserted into the proton self-energy, 
the cut of this amplitude generates five-quark Fock states 
of the proton \uudQQ), see Fig. 

Intrinsic strange, charm, and bottom quarks are thus 
a fundamental property of the wavefunctions of hadronic 
bound states [SHE]- While the extrinsic contributions to 
the heavy quark parton distribution functions (PDFs) 
are most important at low x and depend logarithmi¬ 
cally on the heavy quark mass Mq, the intrinsic heavy 
quark contributions are dominant at high x and depend 
on 1/Mq. Because the extrinsic heavy quarks are gen- 
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Figure 1: Five-quark Fock state \uudQQ) of the proton 
and the origin of the intrinsic sea. 


erated by gluon splitting, their PDFs are always softer 
than those of the parent gluon by a factor of (1 — a:). 
In contrast, the high x intrinsic heavy quark contribu¬ 
tions are kinematically dominated by the regime where 
the \uudQQ) state is minimally off shell, corresponding 
to equal rapidities of the constituent quarks. The result¬ 
ing momentum and spin distributions of the intrinsic Q 
and Q can be distinct, e.g., s(x) ^ s(x) since the comov¬ 
ing uudQQ quarks are sensitive to the global quantum 
numbers of the proton. 

A finite intrinsic charm contribution to the nucleon 
has been extracted from lattice QCD. An analysis by the 
MILC collaboration |9] yields a probability for the charm 
matrix element {N\cc\N) in the range of 5 — 6%, consis¬ 
tent with a four-loop perturbative QCD calculation m- 

While the first experimental evidence of intrinsic heavy 
quarks came from the EMC measurement of the large x 
charm structure function m, a variety of other charm 
hadron and charmonium measurments are consistent 
with the existence of intrinsic charm. Open charm ob¬ 
servables in hadroproduction include forward Ac produc¬ 
tion at the ISR |12Q and asymmetries between leading 
and nonleading charm (D mesons which share valence 
quarks with the projectile and D mesons which do not, 
respectively) measured as functions of xp and px in fixed- 
target experiments, WA89 and WA82 at CERN; E791 
and SELEX at Fermilab, see Refs. UMa and references 
therein. Previous fixed-target J/ij) measurements also 
give indications of important intrinsic charm contribu¬ 
tions, particularly from the nuclear mass, or A, depen¬ 
dence, as measured by NAS at CERN as well as E772 
and, later, E866 at Fermilab, see e.g. US]. Indeed, the A 
dependence, proportional to A“, is quite different than 


^ Similarly, the coalescence of comoving b, u and d quarks from the 
\uudbb > intrinsic bottom Fock state in the proton can explain 
the high xp production of the Ai,{udb) baryon, as observed at 
the ISR [I2]. 


the a ~ 1 expected from extrinsic-type production HZ]. 
At large xp, there are indications of a A?^^ dependence, 
consistent with a nuclear surface-type interaction instead 
of the volume dependence of pQCD. In addition, the NAS 
collaboration measured double J/ij: production at for¬ 
ward Xp in ttA interactions, difficult to explain without 
an intrinsic charm mechanism |18| . All of these observ¬ 
ables can be studied with higher energies and luminosi¬ 
ties at AFTER@LHC, making precision measurements 
possible for the first time. 

In addition to the typical observables for intrinsic 
heavy quarks, these intrinsic heavy quarks also con¬ 
tribute to a number of more exotic observables and inclu¬ 
sive and diffractive Higgs production pp —>■ ppH, in which 
the Higgs boson carries a significant fraction of the pro¬ 
jectile proton momentum [iniiini. There are also impor¬ 
tant implications for intrinsic charm and bottom quarks 
in Standard Model physics, as in the weak decays of the 
H-meson |2T] and a novel solution to the J/ip ^ pir prob¬ 
lem HI]. AFTER@LHC could also shed light on these 
topics. 

The rest of this paper is organized as follows. In Sec.jlT] 
we give an overview of the theoretical models predicting 
the cc-shape (but not the normalization) of the intrinsic 
charm and bottom parton distribution functions. In Sec. 
|HI[ we discuss the constraints on the normalization of 
the intrinsic charm (IC) obtained in global analyses of 
PDFs. Section IV is devoted to the intrinsic bottom (IB) 
content of the nucleon for which there are currently no 
quantitative constraints. In Sec. j^we review collider ob¬ 
servables sensitive to an intrinsic charm or bottom PDF. 
Finally, in Sec. |VI| we present our conclusions. 


II. THEORETICAL MODELS 

The QCD wavefunction of a hadron can be represented 
as a superposition of quark and gluon Fock states. For 
example, at fixed light-front time, a hadron wavefunction 
can be expanded as a sum over the complete basis of free 
quark and gluon states: 1^'^,) = 1"^) ^t,*) 

where the color-singlet states, |m), represent the fluctu¬ 
ations in the hadron wavefunction with the Fock com¬ 
ponents \q 1 q 2 q 3 ), \q 1 q 2 q 39 ), \q 1 q 2 q 3 cc), etc. The boost- 
invariant light-front wavefunctions, ipryifi^{xi,kT^i) are 
functions of the relative momentum coordinates Xi = 
kf /P"*" and kp^i where ki denotes the parton momenta 
and P the hadron momentum. Momentum conservation 
demands = 1 and X]i=i ^T,i = 0 where n is 

the number of partons in state |m). For example, as pre¬ 
dicted by Brodsky and collaborators, in the BHPS model 
intrinsic charm fluctuations 13113 ] can be liberated by a 
soft interaction which breaks the coherence of the Fock 
state m] provided the system is probed during the char¬ 
acteristic time that such fluctuations exist. 

Microscopically, the intrinsic heavy quark Fock compo¬ 
nent in the proton wavefunction, \uudcc}, is generated by 
virtual interactions such as gg -A- QQ where the gluons 
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couple to two or more valence quarks. The probability for 
cc fluctuations to exist in a hadron is higher twist since it 
scales as l/m^ relative to the extrinsic, EC, leading-twist 
production by photon-glnon fnsion |18| . 

The dominant Fock state configurations are not far 
off shell and thns have minimal invariant mass, = 
where m| = mf + {k'^ ■) is the square of the 
average transverse mass of parton i. The general form of 
the Fock state wavefunction for a hadron with mass ruh 
appropriate to any frame at fixed light-front time is 




T{xi,kji_i) 

m\ — 


( 1 ) 


where T is a vertex function, expected to be a slowly- 
varying, decreasing fnnction of m\ — The particle 
distribntions are then controlled by the light-front energy 
denominator and phase space. This form for the higher 
Fock components is applicable to an arbitrary number of 
light and heavy partons. Intrinsic cc Fock components 
with minimum invariant mass correspond to conflgnra- 
tions with equal rapidity constituents. Thns, unlike ex¬ 
trinsic heavy quarks generated from a single parton, in¬ 
trinsic heavy quarks carry a larger fraction of the parent 
momentum than the light qnarks in the state ElUSl- 
The parton distributions reflect the nnderlying shape 
of the Fock state wavefnnction. Assuming it is snfficient 
to nse (fcy) for the transverse momentum, the probability 
distribution as a fnnction of a: in a general n-particle 
intrinsic cc Fock state is 

dPiC ^ ^ 5{l-YTi=iXi) , . 

dx,---dxr. ” ’ 

where Nn normalizes the n-particle Fock state probabil¬ 
ity. 

At LO in the heavy quark limit, fhc, fhc ^ ttih, fhq. 


where 


Q— 

r / \ ^DNAa 

/d(») = T^» 


dt 


[-t -k (mA^ - TOjv)^] „2 


[t - mA] 


F\t) , 
( 6 ) 


with F(t) a form factor at the DNA vertex and = 
— m\ yj(l — y). In this case they chose a monopole 
form factor with A^ = 1-2 GeV. The coupling constant 
was assnmed to be “ —3.795. From heavy quark 

effective theories |13], the c distribution in the D is ex¬ 
pected to be hard becanse in the bound state, the c ex¬ 
changes momenta mnch less than me- They make the 
extreme assnmption that the entire D momentnm is car¬ 
ried by the charm quark, Cjj = x6(x — y). 

Next, Steffens et al. investigated all the charm strnc- 
ture function data with two variants of intrinsic charm 
|27| . The first was that of Eq. Q, called ICl in their pa¬ 
per, while the second was a meson clond model, IC2. In 
the second approach, the c distribution is obtained from 
the light-front distribution of D mesons in the nucleon. 


cl“(a;) 


foi^) 


1 

IbTT^ 



9^ix,kl) 

x)isDA, - 


X 


k"^ + (niA^ — (1 — x)mN)^ 

1 — X 


(7) 


A hard charm momentum distribution was assnmed in 
the D, similar to that of Ref. |25]. The vertex function 
g^{x,k\) is parameterized as + -k 

SjjA ) where is the sqnare of the center of mass 

energy of the DAc system and g^ the conpling constant 
at Sjjj^ = m\. For an intrinsic charm probability of 1%, 
A « 2.2 GeV. The charm distribution is then 


TC2 


(x) 




( 8 ) 


= TV, 6(l 

dxi ■ ■ ■ dXn ” {Xc + XcY ^ 


E-.)' 


(3) 


leading to 


pIC LO 

^2 c 


(x) 


-xc(x) 
9 ^ ' 


9"". 


dx\ ■■■ dXf 


dPic 


dXi ■ ■ ■ dXrdXr 


(4) 


There are many applications of intrinsic charm in 
charm hadron production. See, e.g.. Refs. HMSl EE] 
for more details. 

Paiva et al. have also calculated an intrinsic charm 
component of the nucleon sea within the context of the 
meson cloud model |25| . They assumed that the nucleon 
can fluctnate into DAc- The c distribution in the nucleon 
is then 


xcn{x) 


dy Id (y) 



(5) 


where the charm distribution in the Ac is assnmed to be 
CA, ^ 6{x- 2/3) and /a,(x) = /^(l - x). 

Pnmplin [2S] considered a model where a point scalar 
particle of mass uiq couples with strength g to N scalar 
particles with mass m-i, m 2 , ■ ■ ■, mM- The probability 
density is then 


N 


dP = 


(167r2)^-i(Ar- 2)! 






N 





(s - Sq)^ 
(s - ml) 


\ns)\\ 


(9) 


where sq = ^i)■ The form factor F{s) snp- 

presses higher mass state contributions. If the quark 
transverse momenta are neglected, with me mnch greater 
than all other mass scales, and F{s) = 1, then the BHPS 
model is recovered. Two types of form factors were stnd- 
ied, an exponential |F(s)p = exp[—(s — mQ)/A^], and a 
power law, |F(s)p = l/(s -I- A^)”) where the cutoff A is 
varied between 2 and 10 GeV. 
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Hobbs et al. employed a meson cloud type approach 
but specified the spin and parity of all lowest mass charm 
meson-baryon combinations from the 5-particle \uudcc) 
Fock states of the proton |^. They pointed out that 
treating quarks as scalar point-like particles, as in e.g. 
Ref. |28], does not conserve spin and parity. They cal¬ 
culated the appropriate meson-baryon splitting functions 
for the meson-baryon combinations and found that the 
production of charm mesons would be almost entirely 
through D* mesons. To study the phenomenological dis¬ 
tributions of charm mesons and baryons in this approach, 
they studied exponential and confining vertex functions, 
(X exp[—(s — m|,)/A^] and (s — m\)) exp[—(s — m|,)/A^] 
respectively. They used these results to compare to the 
Ac distribution from the ISR ED] and the A^/A^ asym¬ 
metry from SELEX [3T]. See Ref. [2S| for details. 


III. GLOBAL ANALYSES OF PDFS WITH 
INTRINSIC CHARM 


In the standard approach employed by almost all global 
analyses of PDFs, the heavy quark distributions are gen¬ 
erated radiatively, according to DGLAP evolution equa¬ 
tions [32H31I, starting with a perturbatively calculable 
boundary condition |35[ E5] at a scale of the order of 
the heavy quark mass. In other words, there are no free 
fit parameters associated to the heavy quark distribution 
and it is entirely related to the gluon distribution func¬ 
tion at the scale of the boundary condition. As a conse¬ 
quence, also the PDF uncertainties for the heavy quark 
and the gluon PDFs are strongly correlated as has been 
discussed in the context of inclusive Higgs production at 
the Tevatron and the LHC m- However, a purely per¬ 
turbative treatment might not be adequate, in particular 
for the charm quark with a mass rric ~ 1.3 GeV which is 
not much bigger than typical hadronic scales but also for 
the bottom quark with a mass mb ~ 4.5 GeV. Indeed, as 
discussed above, light-front models predict a nonpertur- 
bative (’intrinsic’) heavy quark component in the proton 
wave-function Emsj. Motivated by the theoretical pre¬ 
dictions of the BHPS light-front model, analyses of the 
charm distribution in the proton going beyond the com¬ 
mon assumption of purely radiatively generated charm 
date back almost as far as the BHPS predictions them¬ 
selves. For definiteness, in the following we refer to the 
radiatively generated charm by co(x, Q) and to the intrin¬ 
sic charm by Ci{x, Q). The full charm parton distribution 
is then given by the sum c{x,Q) = co(x,Q) + ci(x,Q). 
Strictly speaking, this decomposition is defined at the ini¬ 
tial scale Qo — me of the DGLAP evolution but holds to a 
good approximation at any scale since the intrinsic com¬ 
ponent Cl is governed (to a very good approximation) by 
a standalone non-singlet evolution equation [31]. A sim¬ 
ilar decomposition is understood for the bottom quark 
which will be discussed in Sec. EYl 

The BHPS model of the \uudcc) Fock state predicts a 


simple form for F 2 cix), 

X (10) 

^(1 — a;)(l-I-lOx-I-x^)-I-2a:(l-I-a;) Inx . 
o 


If there is a 1% intrinsic charm contribution to the proton 
PDF, Nb = 36. 

Hoffman and Moore incorporated mass effects and in¬ 
troduced next-to-leading order corrections as well as scale 
evolution [39]. They compared their result to the EMG 
F 2 c data from muon scattering on iron at high x and 
with the intrinsic charm contribution added to the 
leading order calculation of F 2 c by photon-gluon fusion. 

A complete next-to-leading order analysis of both the 
‘extrinsic’ radiatively-generated charm component and 
the intrinsic component was later carried out by Harris et 
al. |7]. The EMG data with 17 = Q'^j2mpX = 53, 95, and 
168 GeV were fit by a sum of the extrinsic and intrinsic 
components |7|. The normalization of the two compo¬ 
nents were left as free parameters, 


F2c{x,g?,ml) = eF^l{x, ,ml) + 5F^'^^{x, g? ,ml), 


( 11 ) 


with the scale fi = ^JmF Q'^. The parameter e, typi¬ 
cally larger than unity, was considered to be an estimate 
of the NNLO contribution to the extrinsic contribution. 
Since a 1% normalization of the IG component was as¬ 
sumed in Eq. (Ill, the fitted value of 6 is the fraction 


of this normalization. Given the quality of the data, no 
statement could be made about the intrinsic charm con¬ 
tent of the proton when P = 53 and 95 GeV. However, 
with D = 168 GeV an intrinsic charm contribution of 
(0.86 ± 0.60)% was indicated. These results were consis¬ 
tent with those of the original analysis by Hoffman and 
Moore [39] . 

The BHPS light-front model assumes that Ci(x) = 
ci(x). Meson cloud models, introduced later, treat the 5- 
particle Fock state as a combination of (predominantly) 

D A+. In this case, of course, Ci(x) ^ Ci(x) with the 
c quark in the D carrying more momentum than the c 
quark in the charm baryon. An analysis by Steffens et 
al. in the context of the meson cloud model and using a 
hybrid scheme to interpolate between massless evolution 
at high and ‘extrinsic’ production at low found a 
limit of ^ 0.4% |27| . 

Regardless of whether or not the models predict Ci (x) — 
ci(x) > 0, intrinsic charm should provide the dominant 
contribution to the charm density in the proton at large 


For some time, no other analyses of the charm struc¬ 
ture function were made. The EMG data remain the 
only measurement of the charm structure function in the 
relevant {x,Q‘^) regime and are the only DIS data cited 
as evidence for intrinsic charm. The HERA data on F 2 c 
were at too low x to address the issue. 
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The first global analyses of the proton PDFs with an 
intrinsic charm contribution included were performed by 
members of the CTEQ collaboration iiniiii]. In addi¬ 
tion to the BHPS and meson cloud approaches, they also 
allowed for a ‘sea-like’ contribution with the same shape 
as the radiatively-generated charm distribution. They 
characterized the magnitude of the intrinsic charm com¬ 
ponent (ci(a;, Q^)) by the first moment of the charm dis¬ 
tribution at the input scale Qo = rric = 1-3 GeV(3 

ci{N = l,Ql)= [ dx ci{x,Ql)= 0.01, (12) 

Jo 

which translates into a momentum fraction 

(x)ci-i-ci = [ dxx[ci{x,Ql) -|-ci(x,Qq)] = 0.0057 . 

Jo 

(13) 

They found that the global analyses of hard-scattering 
data provided no evidence for or against the existence of 
intrinsic charm up to (x)ci-i-ci = 0.0057, i.e. the qual¬ 
ity of the fit is insensitive to {x)ci+ci in this interval. 
They also found that the allowed range was greatest for 
the sea-like IC, expected since this shape is rather easily 
interchangeable with other sea quark components while 
the other, harder, charm distributions are not |4D]. In 
addition, they concluded that the enhancement due to 
IC relative to analyses without it persisted up to scales 
of ~ 100 GeV and could have an influence on charm- 
initiated processes at the LHC, as is discussed later. The 
CTEQ6.6C proton PDFs were generated as a result of 
this analysis m- 

There are two recent updates to the global analyses, 
reaching different conclusions about the importance of 
intrinsic charm. The first, by Dulat et al. |35], follows the 
previous work in the context of the CTEQ collaboration 
[40l l4T| . The second, by Jimenez-Delgado et al. |43| . 
included more lower energy data than the previous global 
analyses. 

The result of Dulat et al. |35] was based on the 
CTIO NNLO parton densities. Here the strong coupling, 
as{Q^), the evolution equations and the matrix elements 
are calculated at NNLO. Only the inclusive jet data still 
required NLO expressions. Their analysis included DIS 
data from BCDMS, NMC, CDHSW, and CCFR; SIDIS 
data from NuTeV and CCFR; the combined DIS and F 2 c 
data from HERA; Drell-Yan production; the W charge 
asymmetry and rapidity from CDF and DO; and 
the inclusive jet measurements from CDF and DO, see 
Ref. [12] for a complete list. 

Two models of IC were considered: the BHPS light- 
front model and the sea-like IC introduced in Ref. m- 
They found a broader possible probability range for IC 


^ Note that at Qo = the radiatively generated charm com¬ 
ponent {co{x,Q‘^)) vanishes at NLO in the MS scheme so that 
c(a;, Qg) = ci{x,Ql). 


in this analysis, (x)ic = (x)ci+ci(Qo) ^ 0.025 for BHPS 
and (x)ic ^ 0.015 for the sea-like IC, summarized in 
Fig-i This finding differs from the previous work which 
found a larger upper limit on IC for the sea-like model. 
They believe that the difference is caused by the im¬ 
proved treatment of the charm quark mass in the later 
study |12] . 



Figure 2: (Color online) The global chi-square function 
versus charm momentum fraction (x)ic- The two curves 
are determined from fits with many values of (x)ic- 
Two exemplary fits for each IC model are shown as 
dots. Blue dots denotes the BHPS model; the dots have 
(a:)ic = 0.57% and 2%, which are denoted as BHPSl 
and BHPS2. Red denotes SEA model; the dots have 
(j^)ic = 0.57% and 1.5%, which are denoted SEAl and 
SEA2. Additionally the dotted lines show global 
chi-square function with additional penalty, T 2 (i), used 
to set the upper limits on the allowed IC component. 

(Figure taken from |42| i 

In addition to the global fit, they also tested the sen¬ 
sitivity of their result to individual experiments by in¬ 
troducing a penalty factor, T 2 {i), for each experiment i. 
This penalty factor is designed to increase more rapidly 
than the Xi for that experiment when xf goes beyond 
the 90% confidence level. The penalty factor employs 
an equivalent Gaussian variable Sn which measures the 
goodness of fit for each individual data set. Values of 
Sn < |1| are considered good fits, 5'„ > 3 is considered 
to be a poor fit, and values of Sn < —3 are better fits 
than expected from usual statistical analyses. Using the 
Sn dependence on (x)ic, they determined which of the 
data sets used in the global analyses are most sensitive 
to intrinsic charm. The upper limit on the BHPS value 
of (x)ic comes from the CCFR structure function data 
while the HERA combined charm data sets the upper 
limit on IC from the sea-like model m- 

They also studied the sensitivity of their sea-like result 
to the charm quark mass and found that, if the charm 
quark mass was raised from 1.3 GeV, as in the CTIO fits, 
to 1.67 GeV, then the minimum for the global analyses 
would support (x)ic = 0.01 rather than 0 although the 
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global is worse for the larger charm mass |35]. Finally, 
they showed how W and Z production at the LHC might 
be affected by a nonzero IC contribution. 

In the most recent study, Jimenez-Delgado et al. ^5] 
included the full range of high energy scattering data by 
using looser kinematic cuts > 1 GeV^ and W'^ > 3.5 
GeV^. In particular, they included the lower energy 
SLAG fixed-target data which did not pass the more 
stringent standard DIS cuts on the plane ap¬ 

plied in the previous work [301-02] • The EMC F 2 c data, 
cited as the strongest evidence for intrinsic charm in DIS, 
are used as a consistency check. The low energy, high¬ 
er, fixed target data lie precisely in the region where IC 
is expected to be most important. Thus including these 
data could enhance the sensitivity of the global fit to 
IC. Note, however, that some of these newly-added data 
are on heavier targets than the deuteron and thus tar¬ 
get mass corrections, nuclear corrections for A > 2, and 
higher-twist effects need to be taken into account |43| . 

They followed the framework of the JR14 [33] global 
fit which decomposed F 2 into light and heavy compo¬ 
nents. The charm component is itself separated into the 
’extrinsic’ and intrinsic charm components. The fixed- 
flavor number scheme is used to compute the extrinsic 
contribution. In this scheme, the charm quark mass en¬ 
ters the PDF evolution only indirectly through the run¬ 
ning of as |43j . They employed a charm quark mass of 
1.3 GeV, as did Dulat et al. [32] • They used all three in¬ 
trinsic charm models previously considered: BHPS, the 
meson-cloud model (this time including pseudoscalar and 
vector mesons as well as spin 1/2 and spin 3/2 charm 
baryons - the CTEQ analyses only included the scalar 
DAc fluctuation), and the sea-like component |32|- The 
IC contribution was evolved up to NLO. 

They found that the total is minimized for (a;)ic = 
0 with (a;)ic < 0.1% at the 5cr level. When a hadron 
suppression factor to suppress charm contributions near 
threshold is applied, they find a minimum Eit {x)ic = 
(0.15 ±0.09)% for the full data set. The SLAG F 2 (large 
x), NMC cross sections (medium x) and HERA ±20 (low 
x) display the greatest sensitivity to IC, see Fig. for 
details. However, fits without the SLAG data still give 
a low IC contribution [33] . The difference between their 
results and previous results is in part due to the very 
different tolerance criteria, Ax^ = 1 for their fit and 
Ax^ = 100 for Dulat et al. [32] • Increasing the tolerance 
to Ax^ = 100 would also accommodate (a:)ic = 1% at 
the Icr level |43| 

When checked against the EMC F 2 c data, a clear pref¬ 
erence for IC is found, as expected, for the highest-a; data. 
However, these data are typically not included in global 
analyses due to their greater tension with other global 
data sets. 


^ For a critical discussion of the analysis in m and in particular 
of the tolerance criterion Ax^ = 1 see Ref. m- 





Figure 3: (Color online) Contributions to the total x^ 
(black circles), relative to the value Xo for no IC, of 
various data sets as a function of the momentum 
fraction {x)ic- 
(Figure taken from [43 j l 


Given that the two most recent analyses set signifi¬ 
cantly different limits on IC, it is important to collect 
further large-a; data, particularly on F 2 c to try and place 
greater confidence on the limit of IC in the nucleon. 
This would be an important measurement at the future 
electron-ion collider. 


IV. PREDICTIONS FOR INTRINSIC BOTTOM 


In contrast to the case of intrinsic charm, there is cur¬ 
rently no global analysis available that investigates the 
possibility of an intrinsic bottom (IB) content of the nu¬ 
cleon. The main reason for this is the lack of experimen¬ 
tal data that could constrain it. The BHPS light-front 
model |5] predicts the existence of IB with an a;-shape 
very similar to the one of IC given in Eq. (101 but with a 
normalization which is parametrically suppressed by the 
ratio rn^/m\. This fact, together with the observation 
that the IB PDF is governed (to an excellent approxima¬ 
tion) by an independent non-singlet evolution equation 
[55] , can be used to investigate IB in a flexible way with¬ 
out the need of a dedicated global analysis. Such a study 
has been done in Ref. [38] where a set of decoupled IB 
(and IC) PDFs has been provided and used together with 
the CTEQ6.6 PDFs [31] to estimate the impact of the IB 
on new physics searches at the LHC. The advantage of 
this approach is that the provided IB (IC) PDF can be 
used with any standard set of PDFs and the normaliza¬ 
tion of the intrinsic component can be freely adjusted. 
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This is especially useful for studies of possible IB effects, 
as in that case, there are no experimental limits on what 
amount of IB is allowed. 

In the following we show some of the results found 
in Ref. |3S]. In this work, the boundary condition for 
the IB distribution was modeled using the IC distri¬ 
butions in the CTEQ analyses 001 m] scaled down by 
the mass factor The result of such an intrin¬ 

sic bottom distribution 6 i(a;, Q^), with normalization 
/p dxbi{x, ml) = 0.01 is shown in Fig. 4 where 

the ratio of the intrinsic (&i) and the radiatively gener¬ 
ated (bo) component of the bottom PDF is plotted. As 



Figure 4: (Color online) Ratio of intrinsic ( 6 i) and 
dynamically generated (bo) bottom PDFs for various Q 
scales. The perturbative bottom PDF from CTFQ 6 . 6 cO 
[41] is used, the normalization of the IB is taken to be 
such that /p dxbi{x,m1) — 0.01 x 
(Figure taken from |;i 8 [ l 

always in the light-front models the intrinsic component 
is mostly present at large x values. We can see that for 
low scales Q ~ 10 GeV the modification of the bottom 
PDF, Kf, = 1 + bi/bo, can reach = 2.5. However, it 
decreases rapidly with the rising scale. Since bi evolves 
independently of the other PDFs the change in the nor¬ 
malization of the IB component in Fig. can be done by 
simply rescaling the curves in the figure. If we allowed for 
a 0.035 X m^/ml normalization of the IB the modification 
of the bottom PDF would be given by Kf, = I-I- 61/60 x3.5, 
which for high x and Q ^ 10 GeV would result in an en¬ 
hancement of the bottom PDF by a factor ~ 6.25. How¬ 
ever, at a scale of around 100 GeV and x below 0.2-0.3, 
even with the higher IB normalization, the effect is be¬ 
coming negligible. 

In Fig. we show the sum of the intrinsic bottom PDF 
bi and the dynamically generated PDF bo from GTEQ 6.6 
for different normalizations of the IB component, namely 
0.01 and 0.035 xm^/m^. We compare this sum to the 
asymmetric uncertainties of the GTEQ 6.6 PDF set (up- 


The asymmetric errors are computed following Ha HZ]. 


per panel). In the same figure is also shown the ratio of 
the same PDFs to the central value of CTEQ 6.6 (lower 
panel). As can be seen, the IB curve with the 0.035 
xm^/ml normalization clearly lies outside the uncer¬ 
tainty band whereas the one with the smaller normal¬ 
ization is marginally outside the band (up to x 0 . 6 ). 



X 


Figure 5: (Color online) CTEQ 6.6 + bi for different 
normalizations of the intrinsic bottom-quark PDF at 
the scale Q = 10 GeV, compared to the asymmetric 
PDF errors from the same set (upper panel). Also 
shown is the ratio of the same PDF sets to the central 
value of CTEQ 6.6 (lower panel). 

If we are looking for new physics with couplings pro¬ 
portional to the mass, the suppression of IB compared 
to the IC would be partly compensated by the square of 
the coupling. For a more detailed study of the relevant 
parton-parton luminosities please see Ref. |38| . 

V. COLLIDER OBSERVABLES 

Several collider observables receive large contributions 
from heavy quark initiated subprocesses and are hence 
potentially sensitive to an intrinsic charm content in the 
nucleon. In order to expect optimal effects the heavy 
quark PDF should be probed at large x > 0.2 (for light- 
front models) and not too large factorization scales. This 
kinematic region is best accessible at lower energies in 
the center-of-mass system (cms) and/or large rapidities. 
Therefore, a fixed target experiment like AFTER@LHC 
m operating at a cms energy y/s = 115 GeV with a 
high luminosity is ideally suited for searches of IC effects. 
In the following we review some of the collider processes 
which have been studied in the literature in this respect. 

A. Open heavy flavor production 

Inclusive charm hadron (D°, Ac,...) produc¬ 

tion in hadronic collisions was advocated in Ref. 05] as 
a laboratory to probe IC inside the colliding hadrons. In 










































this analysis, predictions for the differential cross section 
in dependence of the transverse momentum pT were ob¬ 
tained in the general-mass variable-ffavor-number scheme 
(GM-VFNS) |491 - l^ at next-to-leading order (NLO). In 
this scheme, the charm quark is an active parton and the 
differential cross sections of inclusive charm meson pro¬ 
duction depend heavily on the PDF of the charm quark. 
The sensitivity of these cross sections to IC was stud¬ 
ied for the Tevatron at a cms energy of 1960 GeV and 
the Relativistic Heavy Ion Gollider (RHIG) at cms en¬ 
ergies of 200 GeV (RHIG200) and 500 GeV (RHIG500). 
The different IG models from the GTEQ6.5c global anal¬ 
ysis uni were employed together with the fragmentation 
functions for charm mesons from Ref. |52]. While the 
effects at the Tevatron were found to be very moderate 
and likely not testable, large enhancements were found 
at RHIG200 reaching values of ~ 3 at pr = 20 GeV. Un¬ 
fortunately, the measurements at RHIG200 are limited 
by the luminosity. At RHIG500 the cross section is in¬ 
creased by about a factor 3.6. However, the sensitivity 
to IG for the light-front models is greatly reduced. 

More recently, the GM-VFNS was applied to obtain 
predictions for the production of inclusive D mesons at 
the LHG for a cms energy of 7 TeV (LHG7) [^. It 
was found that the production cross sections at large ra¬ 
pidities 2 / > 4 are sensitive to an IG component. These 
predictions can be tested by measurements at forward 
rapidities with the LHGb detector. 

The ideal experiment to search for the effects of IG 
would be a high luminosity fixed target experiment such 
as AFTER@LHG operating at a cms energy of 115 GeV. 
In Fig. [^we show results for inclusive D* meson produc¬ 
tion as a function of the transverse momentum of the D* 
meson and integrated over the rapidity range 2 < y < 5 
(in the laboratory frame) in essentially the same setup as 
in Ref. |3H] to which we refer for details. The only dif¬ 
ference is that, following Ref. (54], the default choice for 
the renormalization and factorization scales is = rriT, 
Pf = I^'f = 'mTl2 where tot = \Jpt is the trans¬ 

verse mass. The theoretical predictions are shown on 
an absolute scale in Fig. (left) and as a ratio with re¬ 
spect to the default results in Fig. (right). In both fig¬ 
ures, the black dotted lines have been obtained by vary¬ 
ing the renormalization scale around the central choice 
to pfi = fnTl2 (upper line) and = 2 m 7 ' (lower line). 
In the right figure we repeat the calculation of the cen¬ 
tral prediction in turn with PDF sets GTEQ6.5Gn for 
n = 1 ,... ,6 and normalize the outcome to the default 
prediction with zero IG of Fig. (left). We observe that 
the ratios for n = 1,2, 3,4 corresponding to the BHPS 
(n = 1, 2) or meson-cloud (n = 3,4) models become very 
large at large pt- Indeed, the default cross section can be 
increased by more than a factor 5 at = 20 GeV in sce¬ 
narios with maximally allowed intrinsic charm (n = 2,4). 
Even for the IG sets with smaller normalization (n = 1, 3) 
corresponding to {x)ci+ci = 0.57% and {x)ci+ci = 0.96% 
the cross section would be enhanced by a factor larger 
than 2 (red solid line) or 3 (blue dashed line) at px = 20 


GeV. It is also interesting to note that the phenomenolog¬ 
ical models for a sea like IG (n = 5, 6 ) lead to a significant 
enhancement of the cross section at small px ~ uic which 
would be probed at AFTER@LHG as well. 


B. Production of a photon in association with a 
charm quark 

Another process with a wide range of phenomenolog¬ 
ical applications in pp, pA, and AA collisions [SSHSTj 
which is very sensitive to the heavy quark PDF is the 
associated production of a photon with a heavy quark. 
A dedicated study of this process at the LHG operat¬ 
ing at = 8 TeV (LHG 8 ) was performed in Refs. 
(581 [55] where it was demonstrated that the existence 
of IG in the proton can be visible at large transverse mo¬ 
menta of the photons and heavy quark jets at rapidities 
1.5 < |y.y| < 2.4, Ij/cI < 2.4. Indeed, for the BHPS model 
the cross section can be enhanced by a factor of 2-3 for 
pi^ > 200 GeV (see Fig. 5 in [55]). This comes with the 
penalty that the cross section falls rapidly with increas¬ 
ing transverse momentum so that this measurement will 
be limited by statistics. 

Again, as for open heavy flavor production, the lower 
cms energy together with the high luminosity makes 
a fixed target experiment like AFTER@LHG the ideal 
place to discover IG using 7 -f c production. This can be 
seen in Fig. ]^ where the differential cross section is en¬ 
hanced by a factor 5 at p^ = 20 GeV (right panel) with 
a not too small cross section (left panel). 


C. Vector boson production 


Dulat et al. (4^ studied the sensitivity of and 
production to the presence of IG. Vector boson produc¬ 
tion at the LHG is an interesting ground for IG because 
it is at relatively large x for colliders and —>■ is a 

rather clean final state. They did a NNLO calculation of 
W and Z production including IG based on their global 
fits at ^/s = 8 and 14 TeV. They also studied the ratio 
d<7yF++w-{y)ld(Tzo{y) relative to the result with no IG. 
Neither of these calculations showed an effect larger than 
the uncertainties due to the GTIO sets themselves. How¬ 
ever, when the Z^ px distribution with IG was compared 
to that without, they saw a factor of two enhancement 
at Px ^ 500 GeV for = 8 TeV in the range \r]\ < 2.1. 
The corresponding enhancement at 14 TeV was smaller 
at the same px because the x value reached is reduced at 
the higher energy |42[ . 

We show a simple test case here for W and Z produc¬ 
tion to NLO at vs = 7 TeV. We use only the BHPS IG 
parame teri zation for the five-particle Fock state, shown 
in Eq. (10). We assume a 1% normalization and no 


evolution to maximize the possible effect at forward ra¬ 
pidity. The p^-integrated rapidity distribution is shown 
in Figj^ as is the ratio of the result with IG to that with- 
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Figure 6: (Color online) NLO predictions for inclusive D* meson production at AFTER@LHC vs the transverse 
momentum of the D meson. (Left) Differential cross section on an absolute scale without intrinsic charm. (Right) 
Ratio w.r.t. to the central prediction of the left plot. Shown are results using the IC parametrizations from Ref. m 
for n = 1 (red, solid line), 2 (violet, dotted line), 3 (blue, dashed line), 4 (green, long dashed line), 5 (cyan, 
dot-dashed line), 6 (orange, double-dot-dashed line). In both figures, the black dotted lines have been obtained by 
varying the renormalization scale around the central choice (fin = nriT) to = 7X1x12 (upper line) and = 2mT 

(lower line). 




Figure 7: (Color online) NLO predictions for the production of a prompt photon in association with a charm quark 
jet in pp collisions at AFTER@LHC vs the transverse momentum of the photon. Shown are results for an BHPS 
and a sea like intrinsic charm using the CTEQ6.6c PDFs. For comparison, the predictions without an IC using the 
CTEQ6.6M PDFs are shown as well together with the uncertainty band obtained by varying the central 
factorization scale pp = p^ a factor 2 up and down (blue, dotted curves). The right panel depicts the ratio of the 
curves in the left panel with respect to the central prediction without intrinsic charm. 
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out as a function of rapidity. The rapidity distributions 
without IC are given by the solid curves while the dashed 
curves are the calculations with the BHPS IC contribu¬ 
tion to the charm parton density. With BHPS IC, one 
expects enhancement only at forward rapidity. The en¬ 
hancement from IC appears for \y\ > 2.5. Note that if 
the sea-like IC would be used instead, the enhancement 
would be small but finite over all rapidity. 

The W'^ cross section is largest and most forward 
peaked, because of the ud contribution. The contribu¬ 
tion from the cd part is a very small addition since the u 
valence contribution is large and peaks at large x, mak¬ 
ing the y distribution larger at \y\ ^ 2 than at j/ = 0. 
Indeed, it gives the smallest IC contribution. The W~ 
distribution should have the largest possible contribution 
from IC because both the du and dc peak at low x and 
because the d valence distribution peaks at lower x so 
that the W~ rapidity distribution has a maximum at 
y = Q. At lyl ~ 4, the IC enhancement is ^ 40%. Fi¬ 
nally, the distribution, with a plateau over \y\ < 1.5, 
also has a very small IC contribution because the charm 
enhancement only comes through cc. 

Such IC enhancements are only visible outside the 
midrapidity acceptance of the collider detector coverage 
of CMS and ATLAS. However, LHCb or ALICE cover 
this forward rapidity range with muons and could detect 
forward They could also look at the lepton rapidity 
asymmetry, {kF+ — W~)/-\- W~), at forward rapid¬ 
ity. The statistical accuracy of the measurement would 
need to be high to distinguish an IC enhancement from 
the no IC result, especially since the 1% BHPS IC is 
likely an upper limit on this enhancement. Note that the 
higher energy of LHC Run 2 will reduce the potential 
enhancement even though it would increase the rates. 

VI. CONCLUSIONS 

The existence of non-perturbative intrinsic charm and 
bottom components is a fundamental prediction of QCD. 


In this article, we have reviewed the current status of 
our understanding of this intrinsic heavy quark content 
of the nucleon which yet remains to be confirmed ex¬ 
perimentally. In particular, after introducing theoretical 
models predicting the intrinsic heavy quark distributions 
we have turned to a summary of the available informa¬ 
tion on intrinsic charm coming from global analyses of 
parton distribution functions. There are no global anal¬ 
yses of intrinsic bottom available and we have described 
how IB can be modeled in order to explore its impact on 
collider observables keeping in mind that bottom quark 
initiated subprocesses play an important role in certain 
electroweak observables and in models for physics beyond 
the Standard Model. We then have turned to a discussion 
of collider processes where IC could be discovered. Gen¬ 
erally, the effects of IC are larger at colliders with a lower 
center-of-mass energy and for hard processes with mod¬ 
erate factorization scales. Therefore, a high-luminosity 
fixed target experiment like AFTER@LHC operating at 
a center-of-mass energy = 115 GeV would be ideally 
suited to discover or constrain IC. 
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